Design of a gap tunable flux qubit with FastHenry
Akhtar Naheed, Zheng Yarui, Nazir Mudassar, Wu Yulin, Deng Hui, Zheng Dongning†, , Zhu Xiaobo‡,
Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China

 

† Corresponding author. E-mail: dzheng@iphy.ac.cn

‡ Corresponding author. E-mail: xbzhu16@ustc.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11374344, 11404386, and 91321208), the National Basic Research Program of China (Grant No. 2014CB921401), and the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB07010300).

Abstract
Abstract

In the preparations of superconducting qubits, circuit design is a vital process because the parameters and layout of the circuit not only determine the way we address the qubits, but also strongly affect the qubit coherence properties. One of the most important circuit parameters, which needs to be carefully designed, is the mutual inductance among different parts of a superconducting circuit. In this paper we demonstrate how to design a gap-tunable flux qubit by layout design and inductance extraction using a fast field solver FastHenry. The energy spectrum of the gap-tunable flux qubit shows that the measured parameters are close to the design values.

1. Introduction

Superconducting quantum circuits based on Josephson junctions[1] are promising candidates for the implementation of quantum computers.[2] This is due to the fact that the complex superconducting quantum circuits can be freely constructed by using well-established micro-fabrication technology in the semiconductor industry. Also the superconducting quantum circuits have advantages, such as being easy to control, scalability, and strong inter-coupling as compared with other quantum systems such as ion traps,[3] nuclear spins,[4,5] electron spins,[6] neutral atoms,[7] and photons.[8] The basic component of a quantum computer is the qubit, which is a two-level system. For a superconducting circuit to function as a qubit, its circuit parameters must be well designed. One parameter which can play an important role in affecting the characteristics of a superconducting qubit is the inductance. In this paper, we demonstrate how to choose proper values for the inductances in designing a gap-tunable flux qubit.[9,10]

A superconducting flux qubit loop is interrupted by three Josephson junctions,[11] the critical current of one of the junctions is α times smaller than those of the other two junctions, which are identical and characterized by the Josephson energy EJ and the charging energy Ec. When the magnetic flux threading the qubit loop Φɛ is around half flux quantum Φ0/2, the qubit ground and excited states correspond to clockwise and anticlockwise persistent currents ±Ip in the qubit loop. When the qubit is biased at the symmetric point or the optimal point Φɛ = Φ0/2, the persistent current states are degenerate, quantum tunneling lifts this degeneracy and forms an energy gap Δ. The Hamiltonian of a flux qubit is given by:

Here,

Φ0=h/2e is the flux quantum, while σz and σx are Pauli matrices. A flux qubit in most cases works at the optimal point for good coherence properties. Therefore, it is desirable to make the energy gap Δ tunable. In order to do this, one needs to replace the smaller junction of the flux qubit by a dc superconducting quantum interference device (SQUID, called an α loop). By applying bias current in an α bias line which is inductively coupled with the α loop, the critical current of the dc superconducting quantum interference device can be modulated, thus the value of α can be changed, making Δ tunable. Since the working point of a flux qubit is determined by the magnetic flux threading the qubit loop Φɛ, the qubit working point can be changed by applying bias currents in an mw/flux bias line placed in the vicinity of the qubit loop. Microwave for qubit operations is also applied through the mw/flux bias line. The qubit state readout is performed by coupling the qubit to a dc-SQUID to detect the persistent current in the qubit loop. Coupling between two flux qubits or between a qubit and other systems can also be realized through mutual inductance. Thus we can see that for flux qubits the inductance plays an important role. In the following Fig. 1 we show our gap-tunable flux qubit sample with its schematic circuit diagram.

Fig. 1. SEM image of our gap tunable flux qubit sample (a), and schematic circuit diagram (b), showing the qubit loop (green loop), the SQUID loop (red loop), and the control lines. Crosses represent Josephson junctions.
2. Influence of mutual inductance on qubit coherence

One of the biggest challenges of realizing quantum computers with superconducting qubits is coherence. A qubit in a quantum computer needs to perform a huge number of quantum logic operations to solve a practical problem. Like any other quantum system, a qubit undergoes decoherence due to coupling to noise in its environment. Each noise source with a qubit coupled can be represented by a fluctuation term of some parameter λ in its Hamiltonian δ λ and characterized by its quantum spectral density sλ(ω) = 1/(2π) dtδλ(0)δλ(t)〉 e−iωt.[12]

The effect of decoherence can be separated into two effects, i.e., energy relaxation and dephasing, and characterized by relaxation rate Γ1 and dephasing rate Γϕ respectively,[13]

Here is the qubit frequency, Dλ,⊥ and Dλ,‖ are the transverse and longitudinal term of Dλ = (1/h)H/∂λ, which is the sensitivity of the qubit Hamiltonian to the noise source, and represented by λ. The most prominent noise sources of a gap-tunable flux qubit are fluctuations of the magnetic flux bias in the qubit loop fɛ and magnetic flux in the α -loop fα, the sensitivities of the Hamiltonian to these noises are as follows:

At the optimal point where (ɛ = 0), we have Dfɛ = Dfα = 0, Dfɛ and Dfα are determined only by qubit parameters Ip and Δ, thus in order to minimize decoherence, we can engineer the quantum spectral densities Sfɛ (ω) and Sfα (ω). For a superconducting flux qubit, the main noise sources are the readout and control circuits. As we perform qubit state readout by coupling a dc SQUID to the qubit loop, the circuit noise from the readout circuitry is directly coupled to the current through the SQUID bias current Isqb and the circulating current in the SQUID loop Icir. Because the SQUID loop is placed symmetrically around the qubit loop, only the fluctuation in Icir is coupled to the qubit, the coupling of the fluctuation in Isqb to the qubit loop is zero or negligibly small. The quantum spectral density of the noise from the readout circuit is given by[14]

where Zsq (ω) is the impedance of the readout circuit, Φsq is the magnetic flux threading the SQUID loop, Msq is the mutual inductance between the SQUID loop and the qubit loop. Current fluctuations in the mw/flux bias line and bias line are also coupled to the qubit, and their quantum spectral density is given by[14]

where Zmw (ω) and Zα (ω) are the impedances of the mw/flux bias line and α bias line respectively, Mmw is the mutual inductance between the mw/flux bias line and the qubit loop, Mα is the mutual inductance between the α bias line and the α -loop. We can see that the quantum spectral densities of all the circuit noise sources are proportional to the squared mutual inductances. Thus in order to minimize the qubit decoherence, mutual inductances between the qubit and the control or readout circuits should be as small as possible.

3. Design of a gap-tunable flux qubit

The readout of a flux qubit state utilizes the factor that the two persistent current states |±Ip〉 form a basis for the qubit states. Measuring the persistent current in the qubit loop corresponds to a projective measurement of the qubit in the persistent current state basis. Since the SQUID loop is inductively coupled to the qubit loop, the qubit state influences the critical current of the SQUID. Any change of the persistent current is reflected in the SQUID flux bias Φsq and modulates its critical current Isq. We apply pulses to the readout SQUID and detect its switching. The switching is a stochastic process due to quantum tunneling. The pulse amplitude at which the readout SQUID switching takes place has a distribution with typical spread about 1% of Isq, and the center of the distribution is slightly smaller than Isq. In order to obtain sufficient visibility in the readout signal, we need to have about 1% of change in Isq in a change of qubit state from |−Ip〉 to |+Ip〉, that is 2Msq Ip ≈ 100. For a persistent current Ip = 500 nA, we have Msq ≈ 20 pH. While we want the mutual inductance Msq to be as small as possible for minimizing qubit decoherence, the requirement for the readout to be efficient sets a lower limit for the mutual inductance Msq, which as shown above is typically several tens of pico-Henry.

Just like the mutual inductance between the qubit and the readout SQUID, mutual inductances between the qubit and control lines also have their lower limits. The mw/flux bias line is used for providing fast flux bias to shift the qubit frequency during qubit operations and readout, the maximum amount of this shift is about 30 GHz typically. The maximum current can be applied to the control lines without inducing strong crosstalk or risking the breaking down of superconductivity, and is about 1mA in most samples. Since δ fqb ≈ 2Ip Mɛ Iɛ, to induce a qubit frequency change of 30 GHz for a flux qubit with Ip = 500 nA the mutual inductance between the mw/flux bias line and the qubit loop Mɛ must be at least 0.02 pH. The α bias line is used for tuning the qubit energy gap Δ in a frequency range of several GHz, by numerically calculating the energy gap dependence on the value of α, and we find that this requires a mutual inductance Mα ≈ 0.03 pH between the α bias line and the α loop.

We develop a dedicated layout computer-aided design (CAD) tool to model the qubit circuit, figure 2(a) shows the modeled circuit layout of a gap-tunable flux qubit (qubit a), figure 2(b) shows the modeled circuit layout of a flux qubit coupled to a resonator (qubit b). Black lines represent metal wires with adjustable width and thickness, black dots denote connections or ports. This CAD tool generates a digital layout file for further numerical inductance extraction.

Fig. 2. Circuit layout for inductance extraction. Red crosses denote Josephson junctions. (a) A gap-tunable flux qubit with readout and control circuits (qubit a) and (b) a flux qubit inductively coupled to a resonator (qubit b).

We use the fast field solver FastHenry[15] to extract the inductances numerically. The FastHenry solver is one of the most popular techniques to extract circuit inductances. It uses the fast multipole method (FMM) to accelerate matrix-vector production. The algorithm used in FastHenry is an acceleration of the mesh formulism approach. The linear system resulting from the mesh formulation is solved by using a generalized minimal residual (GMRES) algorithm with a fast multiple algorithm to efficiently compute iteration. The digital layout file is used as the input of FastHenry to extract the inductances. After the calculation, the result is compared with the desired inductance value and the layout is adjusted accordingly. After a few iterations we obtain a layout with inductances close to the desired values. Table 1 shows the inductances of our sample extracted by FastHenry.

Table 1.

Mutual inductances extrated by FastHenry for qubit a and qubit b. Here for qubit b we only extract the mutual inductance between the qubit and the resonator Mqr.

.
4. Experimental results

After the qubit circuit is properly designed the sample is then fabricated and measured. The details about the fabrication and measurement procedures can be found in Refs. [16] and [17]. Figure 3(a) shows the magnetic flux modulation of the readout SQUID critical current Isq. Qubit signals appear as a small shift of the modulation curve δΦsq = 150. Thus the mutual inductance between the qubit and the readout SQUID is Msq = δΦsq/2Ip = 28.4 pH, very close to the designed value. Figure 3(b) shows the qubit energy spectrum, the qubit flux bias Φɛ is applied by the mw/flux bias line. By fitting the spectrum line (red dash line) we obtain Ip = 546 nA. Figure 3(c) shows the Rabi oscillation result measured at the optimal point, the 2-μs decay time of the oscillation, indicating the good coherent property of our sample.

Fig. 3. Measurement results of sample qubit a. (a) Qubit signal. The qubit signal appears as a small shift of the readout SQUID critical current modulation curve, marked by two black dash lines and two arrows. (b) Energy spectrum by applying qubit flux bias with mw/flux bias line. The red dash line is a fit of the spectrum line. (c) Rabi oscillation of the sample, red dots are data.The decay rate of the oscillation is 1.96 μs, extracted by fitting the data with an exponential decay sinusoidal curve (blue curve on the plot), P is the switching probability of the readout SQUID.
5. Conclusions

In this paper, we demonstrate how to design a gap-tunable flux qubit with proper values of mutual inductances between different circuit parts. It is valuable for ensuring long coherence time. It provides enough control freedom over the qubit state and has sufficient readout efficiency. The method demonstrated here can be generalized to other types of superconducting circuits.

Reference
1Josephson B D 1962 Phys. Lett. 1 251
2Nielsen M AChuang I L2000Quantum Computation and Quantum InformationCambridgeCambridge University Press2122–12
3Monroe CKim J 2013 Science 339 1164
4Gershenfeld N AChuang I L 1997 Science 275 350
5Cory D GFahmy A FHavel T F 1997 Proceedings of the National Academy of Sciences of the United States of America 94 1634
6Loss DDiVincenzo D P 1998 Phys. Rev. 57 120
7Deutsch I HBrennen G KJessen P S 2000 Fortschritte Der Physik-Progress of Physics 48 925
8Knill ELaflamme RMilburn G J 2001 Nature 409 46
9Paauw F GFedorov AHarmans C J P MMooij J E 2009 Phys. Rev. Lett. 102 090501
10Zhu XKemp ASaito SSemba K 2010 Appl. Phys. Lett. 97 102503
11van der Wall C Hter Haar A C JWilhelm F KSchouten R NHarmans C J P MOrlando T PLloyd SMooij J E 2000 Science 290 773
12Ithier GCollin EJoyez PMeeson P JVion DEsteve DChiarello FShnirman AMakhlin YSchriefl JSchön G 2005 Phys. Rev. 72 134519
13Paauw F G2009“Superconducting Flux Qubits: Quantum Chains And Tunable Qubits”Ph. D. ThesisNetherland: Delft University Technology
14van der Wal C H2009“Quantum Superposition of Persistant Josephson Currents”Ph. D. ThesisNetherland: Delft Univ. Technology
15Kamon MTsuk M JWhit J K 1994 IEEE Trans. Microwave Theory and Techniques 42 1750
16Wu Y LDeng HYu H FXue G MTian YLi JChen Y FZhao S PZheng D N 2013 Chin. Phys. 22 060309
17Wu Y LDeng HHuang K QTian YYu H FXue G MJin Y RLi JZhao S PZheng D N 2013 Chin. Phys. 22 090312